Superballistic spreading of wave packets 
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We demonstrate for various systems that the variance of a wave packet M(t) oc t v , can show a 
superballistic increase with 2 < v < 3, for parametrically large time intervals. A model is constructed 
which explains this phenomenon and its predictions are verified numerically for various disordered 
and quasi-periodic systems. 
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The time evolution of wave packets in one dimensional 
(ID) and quasi-lD lattices is described by the time- 
dependent Schrodinger equation 



.dc^ = g H 

m=n—b 



dt 



nm^m 



(i) 



where c n {t) is the probability amplitude for an electron 
to be at site n, b is the number of channels, and H nm 
is a tight-binding Hamiltonian. When translational sym- 
metry is present, the eigenstates of the system are plane 
waves and the variance of wave packets increases quadrat- 
ically in time (" ballistic spreading" ) . On the other hand, 
since the pioneering work of Anderson [[l], it is known 
that disorder usually tends to suppress propagation and 
leads to localization. In the one-dimensional case, even a 
small amount of disorder leads to localization of all eigen- 
states |],f|, and therefore asymptotically the spreading 
of a wave packet remains bounded. In higher dimen- 
sions a localization to delocalization transition can occur 
which leads to diffusive or subdiffusive spreading of a 
wave packet ||||. In addition, there are quasi-periodic 
systems which even in ID show fractal energy spectra and 
eigenfunctions leading to a power law spreading [Q that 
is reminiscent of anomalous diffusion in classical systems. 

A global characterization of the dynamical evolution 
of a wave packet is provided by its variance 



M{t) = ^V|c„(i)| 2 cxf 



(2) 



Its time dependence gives a quantitative description of 
the dynamics: v — corresponds to localization, v = 1 
to diffusion, v — 2 to ballistic motion, and v £ (0, 2) to 
anomalous diffusion. It was shown that a ballistic upper 
bound, M[t) < A ■ t 2 , exists for all times with a system 
specific constant A M . Although this statement gives no 
restriction on v for finite time intervals, to our knowledge 
in all studies up to now v < 2 was found for any time. 

In this paper we show that there can be superballistic 
spreading with exponents v £ (2, 3] for parametrically 
large time intervals. Examples with v = 3 can be seen in 
Fig. [l] where the variances of wave packets, while staying 
below the ballistic upper bound, show a cubic growth. 



Each system consists of a perfect lattice with a disordered 
region of finite length. Moreover, extending the length of 
the disordered region increases the time interval of the 
cubic growth. 




FIG. 1. Variance M(t) of a wave packet in a ID disordered 
model of size L = 50, 70, 80, 100 (from left to right) with 
a perfect lattice attached to both ends. The dashed lines 
indicate the time range over which the cubic growth appears 
before the asymptotic ballistic spreading sets in. The inset 
shows the cubic growth for a Band Random Matrix model 
(b= 10, L= 100, T = 0.001). 



We will show that this unexpected behavior of the vari- 
ance is associated with the rate T for the emission of an 
electron from the disordered region into the perfect lat- 
tice. This parameter determines the time scales 



l/V 



and 



1 

*off ~ p 



(3) 



when the superballistic growth starts and ends, respec- 
tively. The different T— dependence of t on and t {{ en- 
sures that this time interval becomes arbitrarily large as 
r decreases. We would like to mention that such interme- 
diate superballistic regimes are well known in other con- 
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texts like in hydrodynamic turbulence (Richardson law) 
or in plasma physics, but are unrelated. 

In order to understand the appearance of the cubic 
growth of the variance in Fig. [I] we consider a simple 
probabilistic model: The disordered part is replaced by a 
point source and anything emitted from it moves with a 
constant velocity v modelling the dynamics of a perfect 
lattice. Initially, all probability is confined to the point 
source and decays with a constant rate T, such that the 
probability at the point source is given by 



P(t) = exp(-Tt). 



(4) 



The variance Mpg(f) of the point-source model is then 
given by 



OO ft 

2 



Mp S (t) = / dxx 2 / dt (-P(t )) 8(x -v(t-t)) , (5) 
Jo Jo 

where —Pit) is the flux emitted from the point source. 
Substituting Eq. (0) for P(t) we get 



M PS (t)=v 2 T dte- Tt {t~t) 2 , 
Jo 



which yields after integration 



(6) 



(7) 



As expected, the variance grows quadratically, Mps(t) ~ 
v t 2 , for asymptotically large times. Expanding the ex- 
ponential term in Eq. (Q) one finds 



(8) 



M mt (t) 



dxx z dt (-P(t )) S(x-L -v(t - * -r)). 



(9) 



This reduces to Eq. (||) of the point-source model for 
times t fulfilling the following three conditions: (i) t 3> r, 
(ii) vt > L, and (iii) M cxt (i) > M int (t). These con- 
ditions set the time scale t on at which the cubic law 
may start, in contrast to t on = for the point-source 
model. Conditions (i) and (ii) give the time scales r 
and L/v, respectively. The internal variance M- lrA {t) is 
bounded by M™ t ax < L 2 . Thus condition (iii) together 
with Mext(i) ~ v 2 T/3t 3 valid for t < t {{ leads to the time 
scale (MX/(" 2 r)) 1/3 . The maximum of these three 
time scales defines the onset of the cubic law 



t on = max-jr, — 



v ' \ v 2 r 



} 



(10) 



In order to see the cubic law over a large time inter- 
val t ff ~ 1/r has to be large implying a small T. For 
sufficiently small V the third time scale in Eq. (M) will 
dominate leading to the scaling presented in Eq. (|3|) for 
v = 3. The ratio t f{/t on scales as T~ 2 / 3 and therefore 
can be made arbitrarily large by decreasing T. 
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Under the condition t < 1/T the cubic term dominates 
all higher orders. Thus we find a cubic increase of the 
variance in the point-source model starting from the time 
t on = up to the time t ff ~ 1/r. At the same time scale 
the crossover to the asymptotic ballistic spreading starts, 
as can be seen from Eq. (]?]). 

An intuitive understanding is based on the fact, that 
the linear decrease of P(t) = 1 — Tt for small times is 
responsible for the cubic growth of the variance [Eq. (ph]. 
During this time the norm of the wave packet outside 
the point source increases linearly. This linear increase 
of the norm combined with the usual quadratic increase 
of the variance due to the ballistic spreading yields the 
cubic growth of the variance. 

A more realistic model should take into account the 
length L of the disordered region and the time r at which 
the norm of the wave packet outside the disordered region 
starts to increase linearly. The total variance is M(t) = 
Mhit(t) + M ext (t) , where Mj nt (t) is the contribution of the 
internal region and M ext (t) originates from the perfect 
lattice, which is given by 
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FIG. 2. (a) Norm of the wave packet inside the disordered 
region of the ID model vs. t for L = 80, 90, 100 showing linear 
decays, (b) Decay rate T vs. L following an exponential law 
(dashed line). 



One can extend the above analysis to other dynamical 
exponents v < 3 by embedding the internal region in a 
lattice showing anomalous diffusion. These lattices are 
characterized by a variance scaling as t M , fi G (0, 2). This 
yields the general expression v = fi + 1 for Eqs. (||J|). 
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Again, the intuitive understanding is that the linear in- 
crease of the norm of the wave packet outside the internal 
region increases the exponent by 1. 
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FIG. 3. Times ton (circles) and t ff (squares) vs. V for 
the ID model (filled symbols) and for the Band Random Ma- 
trix model (open symbols, data for t on and t a are shifted to 
account for different prefactors of the power laws for the two 
models) . 

In the remainder of this paper we will give numerical 
evidence supporting the above analysis. We use perfect 
or quasi-periodic lattices with a finite ID or quasi-lD 
disordered region, although there are many other possi- 
ble settings. Our first example consists of a ID disor- 
dered region of size L attached with semi-infinite perfect 
lattices at both ends. In Eq. (|l|) this corresponds to a 
tridiagonal Hamiltonian (b = 1), with H n ^ n ±\ = 1 and 
H nn = except for a region of size L where H nn is a ran- 
dom number. We fix the disorder strength and we use 
sample sizes L — 50, . . . , 100 such that L > l^, where 
loo is the localization length of the corresponding infinite 
disordered system. In all cases the initial S— like wave 
packet is launched in the middle of the disordered region. 
Fig. [j] shows the variance averaged over 10 disorder real- 
izations. It should be noted, that without averaging we 
get the same qualitative behavior. For small times all 
wave packets spread ballistically until the variance starts 
to saturate when M- mt (t) w l^. Then the cubic increase 
of the variance can be observed before finally the asymp- 
totic ballistic spreading sets in. There is a ballistic upper 
bound for all times in agreement with Ref . || . The range 
over which the cubic law holds increases with the size of 
the disordered sample. This is in agreement with the 
predictions of the point-source model: It is applicable 
since the norm of the wave packet inside the disordered 
region initially decays linearly (Fig. |2|a). From this lin- 
ear decay P(t) = 1 — Tt we determine the rate T which 
decreases exponentially with the sample size L (Fig. |]b), 
as expected for the localized regime ||. Thus extend- 



ing the disordered region decreases T and together with 
Eq. (|^) explains the increase of the superballistic time in- 
terval. The time scales t on and i ff are determined from 
the times when M(t) deviates from the fitted cubic in- 
crease by more than 10%. The scaling of the time scales 
t on and i ff with T are shown in Fig. || confirming Eq. (||). 

We now show that our considerations also apply to 
a Band Random Matrix model that describes quantum 
wires [Q. The Hamiltonian matrix H nrn is real and the 
entries are different from zero in a stripe of width b 
around the diagonal (|n — m\ < b), only. The param- 
eter b defines the hopping range between neighbouring 
sites, or, in the quasi-lD interpretation, the number of 
transverse channels along a thin wire 0. The non-zero 
matrix elements are independent Gaussian random num- 
bers with variance 1 within a region of size L and are 
equal to 1 outside this region. The matrix elements H nm 
that couple the sample to the perfect lattice are random 
numbers with variance T such that r ~ T || . In order to 
study an example where the internal variance Mj nt (t) is 
bounded by the sample size rather than the localization 
length, we choose a sample size L < 1^. Thus the ini- 
tial wave packet spreads diffusively over the disordered 
region before it leaks out to the leads. We find a cubic 
increase of the variance (Fig. [I], inset) and by varying the 
coupling strength T we confirm Eq. @ for t on and i ff 
(Fig- 1). 

Finally, we use a Fibonacci chain model ||] outside 
the disordered region, which allows to vary the exponent 
fi of its anomalous diffusion by changing the potential 
strength V |lCj . Fig. [I| shows the superballistic increase 
of the external variance in nice agreement with the ex- 
pected exponent v = ji + 1. 




FIG. 4. External variance M ex t(t) for a ID disordered 
sample (L — 50) attached to Fibonacci chain models with 
V — 0.1, 0.4, and 0.7 (from top to bottom), showing the 
expected power laws (dashed lines). 
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